1 



Glueball matrix elements on anisotropic lattices 



o 



> 
o 

o 

On 
D 



X 



Kentucky Glueball Collaboration: S.J. Dong, a T. Draper a , W.C. Kuo, b K.F. Liu a , C. Nenkov a *, 
M. Peardon a , J. Sloan a , and B.L. Young c t 

a Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506, USA. 
b Department of Physics, National Sun Yat-Sen University, Kaohsiung 804, Taiwan. 
c Department of Physics and Astronomy, Iowa State University, Ames, IA 50011, USA. 

We describe a lattice calculation of the matrix elements relevant for glueball production in J/tp radiative decays. 
The techniques for such a calculation on anisotropic lattices with an improved action are outlined. We present 
preliminary results showing the efficacy of the computational method. 



1. MOTIVATION 

Glueballs, apart from being one of the most 
exotic states predicted from QCD, are crucial in 
completing the picture of the strong interaction, 
since these mesonic states go beyond the standard 
quark model. Unfortunately there is a consid- 
erable controversy as to their phenomenological 
identification. Even for the low-lying states there 
are a lot of experimental candidates: 



/o(1500) from pp annihilation, 
/j(1710) from radiative J/ip 
decay 



• rr 



)++ 



?;(2225) 



/ 2 (1710) and / 2 (2230) from 
radiative J/ip decay 



Recently, there appears to be strong exper- 
imental evidence that glueballs have been de- 
tected 0]. For years radiative J/ip decays, 



J/1> 



7 + G 



(1) 



have been considered Q the best hunting ground 
for low lying glueballs (i.e. those below 3 GeV). 
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The main objective of this work is to mea- 
sure vacuum-to-glueball transition matrix ele- 
ments which are needed for estimating the partial 
widths of J I ip radiative decays into glueballs || . 
We employ an (^(a^-improved SU(3) gauge ac- 
tion on an anisotropic lattice Q , with the tempo- 
ral lattice spacing much less than the spatial one. 
This has been shown to be an efficient means of 
studying the masses of glueballs, as clearer resolu- 
tion of the decay of correlators is allowed. The use 
of this technique in a matrix clement calculation 
seems highly advantageous as the local current 
operator introduces more UV noise into Monte- 
Carlo measurement. 

2. METHOD 

The matrix elements we compute are of the 
form {G\F^F^\0), where G is the 0++,0-+ or 
2 ++ glueball and F^ u is the gluon field strength 
operator. 

Specifically, we would like to extract the matrix 
elements for the three lightest glueball states that 
can be created by dimension-4 gluonic operators. 
These are; 



• 0++ — Scalar 

(G|Tr£ 2 |0) and (G|TrS 2 |0), 

• ^ — Pseudoscalar 
(G\TtE-B\0), 



(2) 



(3) 
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2++ — Tensor 



(G\TrEiEj\0) and (G\1tBiBj\0). (4) 

Previous studies have shown that smooth inter- 
polating fields give better overlap with the glue- 
ball states. For each lattice irrep of interest, we 
compute the correlation matrix for a large ba- 
sis, {4>i}, of smeared and fuzzed operators. For 
every member of this basis, the correlation with 
the appropriate gluon field-strength operators of 
Eqns. (||-||) is also computed. Variational tech- 
niques are then employed to find the coefficients, 
c,-, needed to make an optimal ground-state cre- 
ation operator, Q(t) = Yli c i'l } i(^)- The rele- 
vant correlators required to extract both the glue- 
ball mass and creation matrix element are then 
formed. These are the smeared-smeared operator 



(0|$(i)$(0)|0) 

t— »oo 



)|$|0} 2 
. (0|$|G)(G|$|0)e- MGt 
and the local-smeared operator 

(0|$0)F 2 (0)|0) - (o|$|o)(o|.f 2 |o) 



(0|$|G)(G|F 2 |0)e- 



-M G t 



(5) 



(6) 



where |G) is the lowest glueball state in the 
irrep. The vacuum subtraction denoted here 
need only be performed for the irrep. A 

three-parameter fit to these correlators then ex- 
tracts the glueball mass Mq and the desired glue- 
ball creation matrix element (G|F 2 |0), as well 
as the glueball interpolation field matrix element 
(0|$|G). 

The construction of the appropriate improved 
gluonic operators on the lattice proceeds as fol- 
lows. For a set of prototype small Wilson loop 
shapes (for example, the plaquette, 2x1 rect- 
angle, etc.) a linear combination of the 48 ori- 
entations Pj (some will be degenerate for loops 
with symmetry) of each prototype, P, is taken 
which transforms under the appropriate irrep, R, 
of the cubic group (including parity). The irreps 
we need to include are thus the A±, A{ ,E + and 
T% . The real part of the trace on these loops 
is then taken, as all the glueball states of inter- 
est are positive under charge conjugation. The 



10" 



a 

10" 2 =- □ 



Scalar Glueball 



Ci _ □ 



$ § r- n 



10"° — 



1% ffl 



- <> 



° 1 r\- i 



10" 



10" 



D 



5 10 15 

Time — slice 



20 



Figure 1. (scalar) correlators. The (3>|$), 

(TrS 2 |$) and (TrE 2 \$) data are indicated by □, 
o and O respectively. The (3>|<&) correlator has 
been rescaled by 2 x 10~ 2 for presentation. 



lattice operators are 

C\ R) {P)= £ a<fReTr(l-[/(P J )), (7) 

3=1. ..48 

where the index i runs over the dimension of the 
irrep. Then, using the standard mapping 



J X 



(8) 



the operators, c[ R \p), are expressed at tree- 
level in perturbation theory in terms of the di- 
mension 4 and 6 continuum gluon field operators 
transforming under the same irrep of the cubic 
group. Once this expansion is performed for all 
prototypes, a linear combination of the lattice 
operators that removes the C(a 2 ) discretization 
errors is found. For the scalar and E ++ polar- 
izations of the tensor currents, this search is un- 
necessary, as conveniently, improving the action 
generates the required operators (Tr Ef and Tr 
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Figure 2. A^f + (pseudoscalar) correlators. The 
($|$) and (TrE • _B|$) data are indicated by □ 
and o. 



Figure 3. E ++ (tensor) correlators. The ($|$), 
(TrB 4 Bj|$) and {TtEiEj\<S>) data are indicated 
by □, o and O respectively. The ($|$) channel 
has been rescaled by 10~ 2 . 



3. PRELIMINARY RESULTS 

In Figs. [j], H and || we present results from a 
feasibility study on an 8 3 x 40 lattice at (3 = 2.4 
and anisotropy a s /a t = 5. These run param- 
eters correspond D to a spatial lattice spacing 
of 0.221(11) fm and temporal cut-off of 4.5(2) 
GeV. For all channels considered, a signal for the 
smeared-local correlator can be observed over a 
sufficiently wide range of time slices to extract 
the relevant observables mentioned in section ^. 
Correlators with chromoclectric fields are seen to 
be significantly more noisy, however. 

4. CONCLUSIONS 

A calculation of the matrix elements relevant 
for predicting the branching ratio for radiative de- 
cays of J I if) into glueballs has been proposed as a 
useful handle on the phenomenological identifica- 
tion of the lighter glueball candidates. This decay 
should be a glueball-rich channel, with fairly clear 
experimental signature. We have performed ex- 
ploratory studies of the ++ , h and 2 ++ (in 



the E irrep) correlator channels, on a lattice with 
5:1 anisotropy and with spatial lattice spacing of 
a s 0.2 fm. They show that anisotropic lattice 
techniques will allow a Monte-Carlo calculation of 
the lattice matrix elements with reasonable com- 
putational resources. 

In the future, we will calculate with various 
lattice sizes and spacings and take the infinite- 
volume and continuum limits. We will include 
heavier glueball states and all 20 irreps of 0£ c . 
We will calculate the renormalization constants. 
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